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Abstract 

We study the spectral densities of (pseudo) scalar currents at finite temperature in 
^ ■ general case when mass of two quarks are different. Such spectral densities are necessary 

I for the phenomenological investigation of hadronic parameters. We use quark propagator 

at finite temperature and show that an additional branch cut arises in spectral density, 
• which corresponds to particle absorption from the medium. The obtained results at T — )• 

. limit are in good agreement with the vacuum results. 

PACS: ll.lO.Wx, ll.55.Hx, 12.38.Bx 

_ 1 Introduction 

X' 

' Heavy ion collision experiments provide an opportunity to investigate particle properties in 
the medium. Inspired by these experiments, there is an increasing interest to investigate the 
properties of hadronic matter under extreme condition [I],[2]- In general, the media created 
by these collisions consist different mesons and baryons. Investigation of hadronic properties 
at finite temperature and density directly using the fundamental thermal QCD Lagrangian is 
highly desirable. However, such interactions occur in a region very far from the perturbative 
regime, where the quark-gluon coupling constant becomes large and the perturbative methods is 
not suitable for calculation of these properties. Therefore, we need nonperturbative approaches. 

Some nonperturbative approaches are Lattice QCD, Heavy Quark Effective Theory (HQET), 
QCD Sum Rules, etc. Among these approaches, the QCD Sum Rules method [3j and its ex- 
tension to finite temperature |1] has been widely used as an efficient and applicable tool to 
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investigate the hadronic properties The main elements in QCD sum rules approach are 
correlation function and dispersion relation. The correlation function is an object with dual 
nature. At large negative g^, it can be evaluated by perturbative methods, whereas at positive 
g^, it must be represented in terms of hadronic observables. Dispersion relation allows us to 
link the thermal correlator for positive values of to the one for negative values. Spectral 
functions in different cases were studied in the literature [6]-[T3]. 

In the present work, we investigate the two-point thermal correlator using the real time 
formulation of the thermal field theory [14J. We calculate the spectral densities of (pseudo) scalar 
currents at finite temperature, which are necessary for the phenomenological investigation of 
(pseudo) scalar mesons parameters. We use quark propagator at finite temperature and show 
how an additional branch cut corresponding to particle absorption from the medium arises in 
spectral density. We also compare our results with predictions obtained at zero temperature. 

2 Thermal Spectral Densities of (Pseudo) Scalar Cur- 
rents 

We begin by considering the thermal correlation function. 



where J{x) =: qi{x)Tq2{x) : is the interpolating current that carries the quantum numbers of 
the state concerned and T indicates the time ordered product. Here F = / or ^75 for scalar and 
pseudoscalar particles, respectively. The thermal average of the operator, A = T{J{x)J~^{0)), 
appearing in the above thermal correlator is expressed as 



where H is the QCD Hamiltonian, and /3 = 1/T stands for the inverse of the temperature. 
Traces are carried out over any complete set of states. In the real time version, thermal 
correlator has the form of a 2 x 2 matrix. However, this matrix may be diagonalized, when 
it is expressed by a single analytic function, which determines completely the dynamics of the 
corresponding two-point function [15]. As this function is simply related to any one, say the 11- 
component of the matrix, we need to calculate only this component of the correlation function. 
The 11-component of the thermal quark propagator is a sum of vacuum quark propagator 
expression and a term depending on the Fermi distribution function. 




(1) 



{A) = Tre-^^A/Tre-^^, 



(2) 




(3) 
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where n{x) is the Fermi distribution function, n{x) = [exp(/3x) + Now, we proceed to 

obtain the temperature-dependent dispersion relation. The time ordering product in Eq. (1) 
can be expressed as 

(T(J(x) J+(x'))) = e{xo - Xo){J{x)J+{x')) + e{xo - xo){J^{x')J{x)), (4) 

where 6{x) is step function. 

Using Kubo-Martin-Schwinger relation, {J{xo)J~^{xq)) = {J^{x'q)J{xq + i^)) for thermal 
expectation and making Fourier and some other transformations, we get the following expression 
for the thermal correlation function in momentum space [II] : 

n(|q|.*) = ^ r HMHM{ ^ - ^^HtM). (5) 

where 

M(|q|,go) = / rf^a;e*^-^(J(x)J+(0)). (6) 

In the above transformations, the following standard integral representation for the 9- step 
function is used: 

^^"^^ -'''>-2v^ loo k^^e • 

The imaginary part of the correlation function can be simply evaluated using the formula 

-±- = TiS{x) + zP(i), which leads to [IE]: 

lo s + Qq 



n(g,T)= I ds^^, (8) 



where p{q,T) = ^ImIl{q,T) ta.nh.^ and Ql = —q^. In some cases, the correlation function 
has ultraviolet divergent. If the spectral density does not vanish at s — )■ oo, the dispersion 
integral in Eq. (8) diverges. A standard way to overcome this problem is to subtract first 
few terms of its Taylor expansion at = from n(g, T). The thermal correlation function in 
momentum space can be written as 

n(g, T) = -tj 0^Jr{TS{k)TS{k - g)), (9) 

where F = / and ^75 for scalar and pseudoscalar particles, respectively. Inserting propagators 
from Eq. (3) in Eq. (9) and carrying out the ko integration, we obtain the imaginary part of 
n(g, T) in the following form: 

f (ik 1 r 

Imll{q, T) = -Nc / —7 (w^ - k^ + k ■ q - Uiqo ± 17111712) 

J 87r"= UJ1UJ2 

X [(1 -111-112 + 2nin2)6{qo - Ui - UJ2) - (^1 + ■n2 - 2nin2)6{qo - uJi + CO2)] 
+ {ujf - k^ + k • q + uiqo ± 17111712) 

X [(1 -111-112 + 2^1^2)5(50 + + ^2) - (ni +112- 2^1^2)5(50 + 0J1- UJ2)] , (10) 



where mi and 777-2 are quark masses, ooi — + , 002 — (k-q)^ + m| , rii — n{oJi), n2 — 
n{u!2) and the plus and minus signs in front of mi, m2, correspond to scalar and pseudoscalar 
particles, respectively. The term, which does not include the Fermi distribution functions, show 
the vacuum contribution. Terms including the Fermi distributions depict medium contributions. 
The delta-functions in the different terms of Eq. (10) control the regions of non-vanishing 
imaginary parts of n(g, T) , which define the position of the branch cuts. As seen the term 
including S{qo — Ui — U2) gives contribution when = cui + 002 ■ Using Cauchy-Schwarz 
inequality (EILi «?)(E^"=i &f ) > (E^"=i ai&i)' we see that. 



UJ1UJ2 = \/k^ + mi\J (k-q)2 + ml> |k||k-q| + mim2, (U) 

and for qo — oui + 002, we get, 

q^^ml + k^ + ml + (k-q)^ + 2ujiUJ2 > (mi + ma)^ + q^. (12) 

Therefore, wc obtain the first branch cut, > (mi + 1712)'^, which coincides with zero tem- 
perature cut describing the standard threshold for particle decays. This term survives at zero 
temperature and it is called the annihilation term. On the other hand, the term including 
5{qo — uJi + UJ2) gives contribution when qo — oui — 002- Similarly to the above expression, we 
obtain, 

= m? + k^ + m^ + (k-q)2 - 2a;ia;2 < (mi - m2)^ + q^ (13) 

and therefore an additional branch cut arises at finite temperature, q^ < (mi — 7712)^, which 
corresponds to particle absorption from the medium. It is called scattering term and vanishes 
at T = 0. 

In the following, we restrict our calculations with |q| = 0, when there is no angular de- 
pendence. Note that, with |q| = 0, the value of |k|, fixed by the 5-functions in Eq. (10), is 
the magnitude of three momentum of quark or antiquark in the center-of-mass of the quark- 
antiquark system: 

u2 _ (go ~ + ^2)') (go ~ ~ ^2)') f,.^ 

In the |q| = case, as it is seen, the terms including 5{qQ — uii —UI2) and 5{qQ-\-uJi-\-L02) functions 
in Eq. (10), give contributions at the regions, q^ > {1711+1712) and go < —{1711+1712), respectively 
giving the vacuum cuts. Similarly, the terms including 6{qo—uji+uj2) and 6{qo+uJi—uj2) functions 
in Eq. (10), give contributions at in the regions, < go ^ (^1 — ''^2) and —(mi —1712) < go < 0, 
respectively giving the Landau cuts. After straightforward calculations, we find the vacuum 
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part of the ImIl{q,T) as: 

/mn(go,T = 0) = ^J{qi -m\- mlf - AmimUql - (mi - m^f). (15) 
oTr^Q ^ ' 

Taking into account both branch cuts after some transformations, the annihilation and scat- 
tering parts of spectral density is found as: 

/ / nT ixfs/ m? — m^w f\fs/ m? — m^wi , , 

Pa,pert(s,T) = po(s)[l " n(^(l + ^ ^ ^ )) - n(^(l - ^— ^))], (16) 

for (mi + m2)^ < s < oo, 

Ps.,ert{s,T) = po(s)[n(^(l + ' ^ ' )) - n( - ^(l - ' ))], (17) 

for < s < (mi — 1712)'^, with mi > m2. Here, po{s) is the spectral density in the lowest order 
of perturbation theory at zero temperature and it is given by 

Po{s) = ^q'{s)v-{s), (18) 

where q{s) = s — (mi — m2)^ and v{s) = (^1 — Amim2/q{s)^ ^ . Here n = 3 and n = 1 for 
scalar and pseudoscalar particles, respectively. As it is seen, at T — )■ limit these expressions 
are in good consistency with the vacuum expressions. Moreover, the obtained results are well 
consistent with the existing results in mi = m2 case [I7]-[19] for the scalar and pseudoscalar 
particles. 

As an example, we present the dependence of the annihilation and scattering parts of the 
spectral density for K"^ and particles in Figs. 1 and 2. In numerical analysis, we use the 
values TUs = 0, 13 GeV and rric = 1,46 GeV for the quark masses. As it is clear, in the region 
of the standard threshold for particle decays, the po{s) is replaced by the annihilation term. In 
the case of light mesons, the values of Pa,pert{s,T) considerably differ from those of the po{s). 
However, in the case of heavy mesons, the Pa,pert{s,T) and po{s) values are very close to each 
other. From Fig. 1, we also see that the in light K"^ cases, the medium contributions play 
important role and consist higher percentage of the total value. 

Our concluding result is that the thermal contributions contribute significantly to the spec- 
tral function. 
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Figure 1: The dependence of the spectral density of meson at temperature T — 120 MeV 
on the \fs parameter. 
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Figure 2: The dependence of the spectral density of meson at temperature T = 120 MeV 
on the ^/s parameter. 
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